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Abstract 

We present studies of single-spin asymmetries for neutral pion electroproduc- 
tion in semi-inclusive deep-inelastic scattering of 5.776 GeV polarized electrons from 
an unpolarized hydrogen target, using the CEBAF Large Acceptance Spectrometer 
(CLAS) at the Thomas Jefferson National Accelerator Facility. A substantial sinc/^ 
amplitude has been measured in the distribution of the cross section asymmetry as 
a function of the azimuthal angle (fih of the produced neutral pion. The dependence 
of this amplitude on Bjorken x and on the pion transverse momentum is extracted 
with significantly higher precision than previous data and is compared to model 

calculations. 

Keywords: 

PACS: 13.60.-r, 13.87.Fh, 13.88.+e, 14.20.Dh, 24.85.+p 

In recent years it has become clear that understanding the orbital motion of 
partons is crucial for achieving a more complete picture of the nucleon in terms of 
elementary quarks and gluons. Parton distribution functions have been generalized 
to contain information not only on the longitudinal momentum but also on the 
transverse momentum distributions of partons in a fast moving hadron. Intense 
theoretical investigations of Transverse Momentum Dependent (TMD) distributions 
of partons and the first unambiguous experimental signals of TMDs indicate that 
QCD-dynamics inside hadrons is much richer than what can be learned from collinear 
parton distributions. 

TMDs were first suggested to explain the large transverse single-spin asymmetries 
observed in polarized hadron-hadron collisions. Since then, two fundamental mecha- 
nisms involving transverse momentum dependent distributions and/or fragmentation 
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functions have been identified, which lead to single-spin asymmetries (SSAs) in hard 
orocesses: a) internal quark motion as represented by, e.g., the Sivers mechanism 



quarks in a nucleon that 
6], which correlates the 



, |2|, |3|, 14, |5[ , which generates an asymmetric distribution o : 
is transversely polarized and b) the Collins mechanism 4, 
transverse spin of the struck quark with the transverse momentum of the observed 
hadron. The 'Sivers-type' mechanism requires non-zero orbital angular momentum 
of the struck parton together with initial- or final-state interactions via soft-gluon 
exchange 3|, |4|, |5[. This mechanism involves TMD distributions which describe the 
correlations between the transverse motion of the parton and its own transverse spin 
or the spin of the initial- or final-state hadron, thereby providing unprecedented 
information about spin-orbit correlations. 

Semi-inclusive deep-inelastic scattering (SIDIS) has emerged as a powerful tool 
to probe nucleon structure and to provide access to TMDs through measurements 
of spin and azimuthal asymmetries. A rigorous basis for such studies of TMDs in 
SIDIS is provided by TMD factorization in QCD, which has been established in 

nrfn fi 

Refs. [3, l8|, |9[ for leading twistjj single hadron production with transverse momenta 
being much smaller than the hard scattering scale. In this kinematic domain, the 



SIDIS cross section can be expressed in terms of structure functions 



llj which 



are certain convolutions of transverse momentum dependent distribution and frag- 



mentation functions. The ana 
of fragmentation functions [12 



ysis of TMDs thus strongly depends on the knowledge 
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15 



16]. 



Many different observables, which help to pin down various TMD effects, are 
currently available from experiments such as: 1) semi- inclusive deep- in-elastic scat- 



1 each twist increment above leading twist (twist-2) contributes an extra suppression factor of 



4 



tering (HERMES at DESY 
and Jefferson Lab 26|, 
PHENIX and STAR at 
hilation (Belle at KEK 




22j, COMPASS at CERN |23|,|2J,|25| 



arized proton-proton collisions (BRAHMS, 
351 ]) and 3) electron-positron anni- 



This letter reports measurements of single-spin asymmetries in the production 
of neutral pions by longitudinally polarized electrons scattered off unpolarized pro- 
tons. The helicity-dependent part ((Jlu) arises from the anti-symmetric part of the 
hadronic tensor llll: 
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The subscripts LU specify the beam and target polarizations (L stands for longi- 
tudinally polarized and U for unpolarized), a is the fine structure constant and fa 
is the azimuthal angle between the leptonic and the hadronic planes defined ac- 
cording to the Trento convention 



The kinematic variables x, y, and z are 



defined as: x = Q 2 /2(P 1 -q),y = (Pi ■ q)/(Pi • fci), z = (Pi ■ P)/(Pi • q), where 



Q 2 



(ki — k 2 ) 2 is the four-momentum of the virtual photon, k± (fc 2 ) is the 



four-momentum of the incoming (scattered) lepton, Pi and P are the four-momenta 
of the target nucleon and the observed final-state hadron, respectively, A e is the elec- 
tron beam helicity, 7 = 2Mx/Q, M and are the nucleon and hadron masses, Pt 
is the transverse momentum of the detected hadron (with P T = P T /\P T \), and p T 
and k,T are the intrinsic quark transverse momenta in the distribution function (DF) 
and fragmentation function (FF), respectively. In Eq. [2] we use small and capital 
letters for DF and FF, respectively. The ratio e of the longitudinal and transverse 
photon flux is given by: e = ^^"^i^ - The structure function i 7 ^™*' 1 receives 
contributions from the convolution of twist- 2 and twist-3 distribution and fragmenta- 



tion functions, such as the twist- 2 Boer-Mulders DF ([391 140]). the Collins FF H±, 
and the twist-3 DFs e and g L . The Boer- Mulders DF h\ describes the correlation 
between the transverse motion of a quark and its own transverse spin, while g 1 - can 
be interpreted as a higher twist analog of the Sivers function. Both functions repre- 
sent spin-orbit correlations. The functions — = — — tt-H^- and — = — — tt-D\ 
are interaction-dependent parts of the higher-twist FFs G x and E, respectively, in 
which m q is the quark mass. The quantities fi and D x are the usual unpolarized 
twist-2 DF and FF, respectively. 

The structure function F^^ h in Eq. |2]is higher-twist by nature. Thus, related ob- 
servables such as beam-spin asymmetries in single-pion production off an unpolarized 
target can only be accessed at moderate values of Q 2 . Such higher-twist observables 
are a key for understanding long-range quark-gluon dynamics. They have also been 
interpreted in terms of average transverse forces acting on a quark at the instant 
after absorbing the virtual photon [41] . 

Different contributions to the structure function in Eq. [2]have been calculated, re- 
lated to both internal quark motion and the Collins mechanisms. Sizable beam SSAs 

r— 1 

were predicted for pion production |42j with spin-orbit correlations as the dynamical 
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origin. Within this framework, the asymmetry generated at the distribution level is 
given by either the convolution of the T-odd Boer-Mulders DF h± with the twist-3 



FF E 
D 1 



44(- 



431 ]. or the convolution of the twist-3 T-odd DF g 1 - with the unpolarized FF 



In contrast, calculations based on the Collins mechanism, eH±, predict vanishing 



beam SSAs for neutral pions [45l. I46l. 1471] . The surprising characteristic that favored 



and unfavored Collins FFs are roughly equal in magnitude but opposite in sign, as 



indicated by the latest measurements from HERMES [22], COMPASS [23j and Belle 
{37], put the 7T° in a unique position in SSA studies since the 7r° FF is the average 
of 7r + and 7r~ FFs. Contributions to the beam SSA related to spin-orbit correlations 
could thus be studied without a significant background from the Collins mechanism. 

Measurements of beam-spin asymmetries in the electroproduction of neutral pi- 
ons in deep-inelastic scattering are presented from the E01-113 CLAS data set using 
a 5.776 GeV electron beam and the CEBAF Large Acceptance Spectrometer (CLAS) 



481 ] at Jefferson Laboratory. Longitudinally polarized electrons were scattered off an 
unpolarized liquid- hydrogen target. The beam polarization was frequently measured 
with a M0ller polarimeter and the beam helicity was flipped every 30 ms to min- 
imize systematic instrumental effects. Scattered electrons were detected in CLAS. 
Electron candidates were selected by a hardware trigger using a coincidence of the 
gas Cherenkov counters and the lead-scintillator electromagnetic calorimeters (EC). 
Neutral pions were identified by calculating the invariant mass of two photons 



detected with the CLAS EC and the Inner Calorimeter (IC) 49(. For events with 
more than two photons, the pair-wise combination of all photons was used. In each 
kinematic bin, tt° events were selected by a Gaussian plus linear polynomial fit to the 
two-photon invariant mass distribution (see Fig{T]). In each <ph bin and for each beam 
helicity, the combinatorial background was subtracted using the linear component of 
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the fit, and it events were selected within the invariant mass region defined by the 
mean of the Gaussian ±3cx, as indicated by the vertical lines in FigHJ 
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Figure 1: (Color online) Invariant mass spectrum of the two photon (77) system M 77 in an arbi- 
trarily chosen x, Pt, z and 0/t-bin, fitted by a Gaussian plus a linear polynomial. Vertical black 
lines indicate ±3cr from the mean. 

Deep-inelastic scattering events were selected by requiring Q 2 > 1 GeV 2 and 
W 2 > 4 GeV 2 , where W is the invariant mass of the hadronic final state. Events with 
missing-mass values for the en system that are smaller than 1.5 GeV (M x (eir°) < 1.5 
GeV) were discarded to exclude contributions from exclusive processes. A minimum 
value for the 7r° transverse momentum, Pj< > 0.05 GeV, ensures that the azimuthal 
angle <ph is well-defined. The total number of selected e7r° coincidences was ~ 3.0 x 10 6 
for the presented z range, 0.4 < z < 0.7, which selects the semi-inclusive region [281 ]. 

The beam-spin asymmetry has been calculated for each kinematic bin 

as: 



where P = 0.794 ± 0.024 is the absolute beam polarization for this data set and N^ 
and N~ are the number of 7r°'s for positive and negative beam helicity, normalized 




(3) 
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Figure 2: (Color online) Examples of fits to the Am asymmetry for 0.4 < z < 0.7, 0.1 < x < 0.2 
and 0.2 GeV< Pt < 0.4 GeV using p sin <j)h (solid line) and poshi(^/(l +picos(ph) (dashed 
line). Both fits yield consistent amplitudes and x 2 per degree of freedom (po = 0.0331 ± 0.0034, 
X 2 /ndf= 1.387 and p a = 0.0329 ± 0.0034, x 2 /ndf= 1.31, respectively). Only statistical error bars 
are shown. 

to the respective integrated charges. The number of 7r°'s is estimated by the integral 
of the histogram in the ±3<r range, minus the integral of the linear component of 
the fit. Asymmetry moments were extracted by fitting the (^-distribution of Alu m 
each x and Pt bin with the theoretically motivated function p o srn 0/i- An example 
of this fit is shown in Fig. [2] for a representative kinematic bin. 

In Fig. [HJ the extracted A s ^f moment is presented as a function of Pt for different 
x ranges. The results are summarized in Table [TJ Systematic uncertainties, repre- 
sented by the bands at the bottom of each panel, include the uncertainties due to 
the background subtraction, the event selection and possible contributions of higher 
harmonics. The first two contributions were estimated as the difference between the 
asymmetry moment extracted from data sets obtained with or without background 
subtraction, and by selecting the 7r° from the combination of all photons in an event 
or from events with exactly two photons. The contribution of higher harmonics was 
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Figure 3: (Color online) Asymmetry moment A s ^y h versus Pt for different x ranges and 0.4 < 
z < 0.7. The error bars correspond to statistical and the bands to systematic uncertainties. An 
additional 3% scaling uncertainty arises from the beam polarization measurement and another 3% 
relative uncertainty from radiative effects which are not included in the band. 

estimated by employing the fit functions posmrfth or p sin 0^/(1 +picos0h). The 
contributions from other harmonics such as sin 20^ or 008 20^ were also tested and 
found to be negligible. All the above contributions were added in quadrature. 

An additional 3% scaling uncertainty due to the beam polarization measurements 
should be added to the above-mentioned systematic uncertainties. Radiative correc- 
tions have not been applied. However they have been estimated to be negligible for 



28 



50j with an overall relative accuracy of 3%. 



the sin0/j modulation 

The A B ^ moment increases with increasing Pt and reaches a maximum at 
Pt ~ 0.4 GeV. There is an indication, within the available uncertainties, that the 
expected decrease of A 8 ^f h at larger Pt could start already at Pt ~ 0.7 GeV. As 
a function of x, A s ^j appears to be flat in all Pt ranges shown in Fig. HI Note, 
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Figure 4: (Color online) Asymmetry moment A s ^j versus x for different Pj* ranges and 0.4 < z < 
0.7. The error bars correspond to statistical and the bands to systematic uncertainties. Comple- 
mentary plot of Fig. [3] 

however, that Q 2 varies with x (see Tabled]). 

The measured beam-spin asymmetry moment for 7r ° ap pears to be comparable 



with the 7r + asymmetry from a former CLAS data set 51] both in magnitude and 
sign, as shown in Fig. El For both data sets the average Pt is about 0.38 GeV. Also 
shown are model calculations of A B ^f h , as indicated in the figure (right-hatched and 
which take only the contribution from Collins-effect eH± into 



left-hatchec 
account 45. 



Dands" 



46 



47 



52], suggesting that contributions from the Collins mechanism 
cannot be the dominant ones. In contrast, preliminary calculations of A 8 ^f h for pions 



53], based on the models from Refs. 
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54] . demonstrate a non-zero contribution 



from g L . Because this DF can be interpreted as the higher- twist analog of the 
Sivers function, it underscores the potential of beam SSAs for studying spin-orbit 
correlations. 
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Beam SSAs for charged and neutral pions were also measured by the HERMES 
collaboration at a higher beam energy of 27.6 GeV [20|. After taking into account the 
kinematic factors in the expression of the beam-helicity-dependent and independent 



terms 



ul) 



1 - y + y 2 /2 

CLAS and HERMES measurements are found to be consistent with each other as 
shown in Figs. [6] and [71 indicating that at energies as low as 4-6 GeV, the behavior of 
beam spin asymmetries is similar to higher energy measurements. For comparison, 
CLAS data in the range 0.4 GeV< Pt < 0.6 GeV are used in Fig. [6] and in the range 
0.1 < x < 0.2 in Fig. [7J because these ranges yield average kinematic values similar 
to HERMES. 

The CLAS data provide significant improvements in the precision of beam SSA 
measurements for the kinematic region where the two data sets overlap, and they 
extend the measurements to the large x region not accessible at HERMES. 

In summary, we have presented measurements of the kinematic dependences of 
the beam-spin asymmetry in semi-inclusive 7r° electroproduction from the E01-113 
CLAS data set. The sin0^ amplitude was extracted as a function of x and transverse 
pion momentum P?, for 0.4 < z < 0.7. The asymmetry moment shows no significant 
x dependence for fixed Pt- Note, however, that Q 2 varies with x (see Tabled]). The 
observed asymmetry moment for ir° suggests that the major contribution to the pion 
beam SSAs originate from spin-orbit correlations. 

The results are compared with published HERMES data 20] . They provide a 



significant improvement in precision and an important input for studies of higher- 
twist effects. Measured beam SSA's are in good agreement, both in magnitude and 



kinematic dependences, with measurements at significantly higher energies 20j, |25|. 
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Figure 5: (Color online) The tt° beam-spin asymmetry moment A s ^j h vs. x compared to that of 
7r + from an earlier CLAS measurement 5l| . Uncertainties are displayed as in Fig. |3l For both data 



sets < Pt >~ 0.38 GeV and 0.4 < z < 0.7. The right-hatched and left-hatched bands are model 
calculations involving solely the contribution from the Collins-effect 47 1. 
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Figure 6: (Color online) Asymmetry moment A s ^y h for n° multiplied by the kinematic factor 
< Q > lf{y) versus x from CLAS and HERMES Q- The 0.4 < P T < 0.6 GeV range of the 
CLAS data is used to compare with HERMES, because this yields average kinematics closest to 
HERMES. 
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Figure 7: (Color online) Asymmetry moment A s ™j for tt° multiplied by the kinematic factor 



< Q > lf{y) versus P T from CLAS and HERMES [20( (the same as in Fig. [6]). The O.K x < 0.2 
range of the CLAS data is used to compare with HERMES, as this yields average kinematics closest 
to HERMES. 
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Tabic 1: The asymmetry moments A s ^/ Ph and their statistical and systematic uncertainties at aver- 
age values of Pt, z, x, Q 2 and y. An additional 3% scaling uncertainty from the beam polarization 
measurement and another 3% relative uncertainty from radiative effects should be added to the 
total uncertainty. 
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